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Abstract 

In this work the collapsing process of a spherically symmetric star, made of dust cloud, in the 
background of dark energy is studied for two different gravity theories separately, i.e., DGP Brane 
gravity and Loop Quantum gravity. Two types of dark energy fluids, namely, Modified Chaplygin 
gas and Generalised Cosmic Chaplygin gas are considered for each model. Graphs are drawn to 
characterize the nature and the probable outcome of gravitational collapse. A comparative study 
is done between the collapsing process in the two different gravity theories. It is found that in case 
of dark matter, there is a great possibility of collapse and consequent formation of Black hole. In 
case of dark energy possibility of collapse is far lesser compared to the other cases, due to the large 
negative pressure of dark energy component. There is an increase in mass of the cloud in case of 
dark matter collapse due to matter accumulation. The mass decreases considerably in case of dark 
energy due to dark energy accretion on the cloud. In case of collapse with a combination of dark 
energy and dark matter, it is found that in the absence of interaction there is a far better possibility 
of formation of black hole in DGP brane model compared to Loop quantum cosmology model. 
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1 Introduction 

The most remarkable and significant discovery in recent past in the field of cosmology is the discovery of the 
fact, that our universe is undergoing an accelerated expansion of late [HI2- This event has overturned a few 
stones in the traditional theories of cosmology. Einstein's equation needed serious revisions in order to account 
for the observed cosmic acceleration. 

Some cosmologists succeeded in affecting proper revisions to the Einstein's equation by introducing changes 
in the right side, i.e., introducing the concept of Dark Energy (DE hereafter) [3] component with negative 
pressure which is responsible for the accelerated expansion of the universe. There have been a lot of research on 
DE till date. Many candidates having the potential to play the role of DE have surfaced. Among them the DE 
associated with a scalar field is often called quintessence. In the cosmological context Chaplygin gas was first 
suggested as an alternative to quintessence. The earliest form of this was known as Pure Chaplygin gas [4j[5]. 
Gradually Generalized Chaplygin gas [5J [5] was formulated and finally Modified Chaplygin Gas (MCG) came 
to the foreground. The equation of state followed by MCG [HI QJJ] is given by 

V = Ap M CG 5 , (1) 

Pmcg 

where p and p are respectively the pressure and energy density and < a < 1, and A, B > 0. MCG appears to be 
consistent with the WMAP 5-year data and henceforth support the unified model with DE and dark matter (DM) 
based on generalized Chaplygin gas [IT] [T2l [T3l HH US] . MCG demonstrates an increasing A behaviour for the 
evolution of the universe [4]. Recent developments regarding MCG can be found in [T6] [17\ [TH] [19]. Recent 
supernovae data also favors the two-fluid cosmological model with Chaplygin gas and matter [20] . 
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In this paper in addition to MCG we will study the gravitational collapse of a star in presence of one more 
fluid known as, Generalized Cosmic Chaplygin Gas(GCCG). Pedro F. Gonzalez-Diaz [21] for the first time 
presented the idea of GCCG. His main motivation was to study DE accretion onto Black Hole(BH), and by 
using GCCG as DE he showed that Big Rip, i.e., singularity at a finite time , is completely out of question, 
which was an essential component of the previous models. The Equation of state for GCCG is given by, 



P ~ Pgccg 



C'+{p G + cc G -Cr (2) 



where C = yq^; — 1. Here A' is a constant which can take on both positive and negative values, and 
> lu > — £, £ being a positive definite constant which can take on values larger than unity. GCCG can explain 
the evolution of universe starting from dust era to ACDM , radiation era, matter dominated quintessence and 
lastly phantom era [22] . 

It is a known that general relativity is incomplete as a theory of gravity. The reason being that the 
gravitational collapse of physical matter produces a BH with a singularity inside, which is a point in space- 
time where the curvature and energy density diverge. At this point all mathematical structures break down 
preventing any further analysis beyond the singularity. Herein lies the motivation for extensive research in 
gravitational collapse. Till date a lot of significant work has been done in gravitational collapse starting from 
the pioneering works of Oppenheimer and Snyder |23j . In this paper our motivation is to study the nature and 
outcome of gravitational collapse of a star made up of DM in the background of DE (of different forms), in two 
types of widely known gravity theories. 

The paper is organized as follows: In section [21 the general formulation of the collapsing process of a cloud 
is given. Section [3] and section [4] deals with the details of collapsing procedure in DGP brane model and Loop 
quantum cosmology (LQC) respectively. Finally the paper ends with some concluding remarks in section [6l 

2 General Formulation of the Collapsing process 

The flat, homogeneous and isotropic FRW model of the universe is described by the line element 

ds 2 = dt 2 + a 2 {t) [dr 2 + r 2 (d9 2 + Sin 2 6d(j) 2 )] (3) 
The energy conservation equation is given by 

Pt+3~{ Pt +Pt) = (4) 

with p T = Pm + Pe and p T = Pm + Pe- 

The interaction Q(t) between DM and DE can be expressed as 

PM + 3-^Pm = Q (5) 

Pe + 3-{p E +Pe) = -Q (6) 
a 

Now, if we consider gravitational collapse of a spherical cloud consists of above DM and DE distribution and is 
bounded by the surface £ : r = then the metric on it can be written as 

ds 2 = dT 2 - R 2 (T){d9 2 + Sin 2 6d(j) 2 } (7) 

Thus on X : T — t and R{T) — r^a(T) where R(r. t) = ra(t) is the geometrical radius of the two spheres t, r = 
constant. Also the total mass of the collapsing cloud is given by 



M{T)=m{r,t)\ = Vad 2 



= l -R{T)R 2 (T) (8) 
s z 



The apparent horizon is defined as 

R, a R,p g aP = 0, i.e., r 2 d 2 = l (9) 
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So if T = Tah be the time when the whole cloud starts to be trapped then 



R 2 {T AH ) ^=ria\T AH ) = l (10) 

As it is usually assumed that the collapsing process starts from regular initial data so initially at t = ti (< Tah)-, 
the cloud is not trapped i.e., 

rlh 2 (U) < 1, (r s d(ti) > -1) (11) 

Thus if equation (fTU)) has any real solution for Tah satisfying (ITT1) then black hole(BH) will form, otherwise the 
collapsing process leads to a naked singularity (NS). So the gravitational collapse and consequently the formation 
of a BH solely depends upon the nature of root obtained from equation (fTU|) . If any real solution for Tah exists 
for equation (fTOf then apparent horizon will be formed and thus a BH. If there is no real solution the collapse 
is destined to result in a NS. 



3 Gravitational Collapse in DGP brane Scenario 

A simple and effective model of brane-gravity is the Dvali-Gabadadze-Porrati (DGP) braneworld model [24l[25l 
|2"5] which models our 4-dimensional world as a FRW brane embedded in a 5-dimensional Minkowski bulk. It 
explains the origin of DE as the gravity on the brane leaking to the bulk at large scale. On the 4-dimensional 
brane the action of gravity is proportional to M 2 whereas in the bulk it is proportional to the corresponding 

M 2 

quantity in 5-dimensions. The model is then characterized by a cross over length scale r c — ^jrfr such that 
gravity is 4-dimensional theory at scales a << r c where matter behaves as pressureless dust, but gravity leaks 
out into the bulk at scales a >> r c and matter approaches the behaviour of a cosmological constant. Moreover 
it has been shown that the standard Friedmann cosmology can be firmly embedded in DGP brane. 
The Friedmann equation in DGP brane model is modified to the equation 




# =1 W^ + — + '— I <'.2) 
and 

2H--)h = -H( Pt + Pt ), (13) 

where H = - is the Hubble parameter, pt and pt are the total cosmic fluid energy density and pressure 

M 2 

respectively and r c — is the cross-over scale which determines the transition from 4D to 5D behaviour and 

e = ±1 (choosing M 2 — 8nG = 1). For e = +1, we have standard DGP(-t-) model which is self accelerating 
model without any form of DE, and effective w is always non-phantom. However for e = —1, we have DGP(— ) 
model which does not self accelerate but requires DE on the brane. 

In this section the role of DM and DE will be shown separately during the collapse. Then we will study the 
effect of the combination of DE and DM in the collapsing process, both in the presence and in the absence of 
interaction. 

3.1 Collapse with Dark Matter 

Here pu 7^ 0, pe = Pe = 0. So from conservation equation ([5]), we have 

PM = % (14) 

where Cq is the constant of integration. Using this relation in equation (jT2"j) we get, for DGP(+) model: i.e., 
for e = +1 

va + r 2 a , / 3 / r\ a 3 

a2 + - log ^ + V« + a 3 J 



2 

3a 



2 



ir + Do (15) 
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Fig. 1 Fig. 2 

Fig 1 : The time derivative of the radius is plotted against time, a = 0.5, r = 10, r c = 100 is 
considered. 

Fig 2 : The mass of the collapsing cloud is plotted against time, a = 0.5, r = 10, r c = 100 is 
considered. 



The corresponding expression for DGP(-) model: i.e., for e 

■ ^ log (a% + y/a 



-1 is, 



2 

3a 



\J a + a 3 



a 

y 



t 

"Zr~c. 



Dn 



(16) 



Where D is the integration constant and a — 



4KCa 



Considering DGP(-) model and using the value of pm 



from equation (|14p and using the relation px = Pm + Pe, the expressions for the time derivative of geometrical 
radius R and mass M(T) of the collapsing object are obtained as follows, 



and 



R(T) 



M{T) 



7 3 / 2 



7 3/2\ 



8ri 



(17) 



(18) 



From the above solutions it is evident that as T — > oo, a — > oo, pm —> 0, R — > 0, M(T) — > 0. So we see that 
there is a tendency of matter density being diminished as time passes, and finally it tends towards zero. The 
time for formation of apparent horizon is given by the real root of the equation, R 2 (Tah) = r^a (Tah) = 1 
as given by equation (|10l) . Thus the corresponding expression for DGP brane is, 



,3/2 



iari 



(19) 



The cloud initially starts untrapped, and gradually as T becomes equal to Tah, the apparent horizon starts 
forming and the cloud begins to be shrouded by the horizon. The result is the formation of a BH singularity. 

Normal collapse of matter influences every collapsing particle to move towards the collapsing centre. So 
for a normal matter it is quite natural to have a negative outward velocity at a particular point. With time 
as the collapse centre gathers more mass, the increasing gravitational attraction pulls the rest particles more 
efficiently causing an increasing negative outward velocity. The same here seems to happen for the case of DM. 
As we do not have any previous idea of the time gradient of DM collapse. So we cannot comment whether this 
is normal or something different. 
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Fig. 3 Fig. 4 

Fig 3 : The time derivative of the radius is plotted against time, r = 10, r c = 1000, C\ = 0.05, X% = 
40, X 2 = 0.04, X 3 = 0.4 is considered. 

Fig 4 : The mass of the collapsing cloud is plotted against time, r = 10, r c = 1000, C\ = 0.05, X\ = 
40, X 2 = 0.04, X 3 = 0.4 is considered. 



3.2 Collapse with Dark Energy in the form of Modified Chaplygin gas 

DE in the form of MCG is considered in this section. Here, om = 0, Pmcg — Apucc — • So the solution 

Pmcg 

for density is given by 



Pmcg 



D 



1+A a 3 ( 1+A )( 1+Q ) 
The expressions for relevant physical quantities are 



R(T) = -rs a 



X 1 



Ci_ 

n x 2 



2X 3 



1 e 

4rj? 2r c 



and 



M(T) = -r|a 3 



Ci_ 

n x 2 



2X-, 



1 : J_ 

4r2 2r c 



(20) 



(21) 



(22) 



where X\ = j^j, X2 = 3 (1 + A) (1 + a), X% = 2 (i+ a ) ■ ^he limiting values of the physical 
parameters are as follows: 



easel 

a ->• : pmcg -> 00, for 1 + A > 0; pmcg 



B 



1 + A 



for 1 + A < 



i?^— ^ ,== , /or 1 + A>0; i?->-oo, for 1 + A < 



M(T) 



r 3 a - 3A C 1+ 



i , for 1 + A>0; M(T) -> 00, /or 1 + A < 



a ->• 00 : pmcg — 



1 + A 



case2 



/or 1 + A>0; p M cc -> 00, /or 1 + A<0 
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Fig 5 : The time derivative of the radius is plotted against time, r = 10, r c = 100000, C = 10, B' 
5, Xi = 0.00005, X 2 = 1.5, X 3 = 0.25 is considered. 

Fig 6 : The mass of the collapsing cloud is plotted against time, r = 10, r c = 100000, C = 10, B 1 
5, Xx = 0.00005, X 2 = 1.5, X 3 = 0.25 is considered. 



R 



-oo, 



for 1 + A > 0; R -> 



— rC 1 a 2 



M(T) 



oo, 



for 1 + A> 0; M(T) -> 



G 



for 1 + A < 
/or 1 + A < 



The limiting value of M(T) for 1 + A > is unphysical for a — > oo. This is because the mass of the collapsing 
cloud has to decrease in presence of MCG due to accretion. The cloud will start untrapped at the instant given 
by the real roots of the following equation and gradually start to be trapped, 



r,X 2 



2X 3 



1 

4r? 



e 

2r~< 



= 1 



(23) 



In Fig. 3 and Fig. 4 we find the time gradient for R, i.e., R and mass against time for DE collapse. Here we 
can see that at first the DE collapse requires an outward negative velocity. However, R has a monotonically 
decreasing magnitude which tends to zero with increasing time. This implies, with time the DE collapse 
decreases, the tendency to collapse. Even the final fate might be the fact that there is no ingoing fluid. The 
change of mass is quite amazing in nature as it sharply indicates that matter is flowing outward from collapsing 
centre. 



3.3 Collapse with Dark Energy in the form of Generalised Cosmic Chaplygin Gas 

as given in 



Here DE in the form of GCCG is considered. So, pu = 0, p = —p G ccG 
equation ([2]). The solution for density is given by 



c + iph+ccG-c) 



PGCCG 



C' + il 



B' 



3(l+Qi)(l+uO 

The expressions for relevant physical quantities are 



, B is positive integration constant. 



(24) 



R(T) = -r s a 



1 



\5 



c 



B> 

a x z 



v: 



2 v.: 



1 e 

Arl 2r c 



(25) 
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and 



M(T) = -rla* 



\5 



C'+ll + 



B' 



1 e 

4r? 2r„ 



(26) 



where X[ = yq^j , = 3 (1 + a) (1 + w), X3 = 2 (i+ a ) ■ The limiting values of the physical parameters are 
follows: 

easel 

1 

a^O: pgccg -> 00, /or 1 + uj > 0; pgccg -> (C + 1) 1+ " , for 1 + u < 

fl4 ^= J B' 2 ( 1 +°)( 1 +-> , /or l+w>0 

V3a 



as 



M(T) -> -!—B'^+<*)(i+^ , /or l + u;>0 
case2 

1 

a -> 00 : Pgccg -> 00, for l+w<0; Pgccg -> (C + 1) 1+Q , /or 1 + cj > 
ra 



M(T) 



3 



4 r 2C"TT^ + 3 + V3e I , /or 1 + w > 



ar 
24rl 



4f2(C" + l)^ +3 + V3e) , /or l+w>0 



The corresponding expression for the formation of event horizon is, 



1 



C' + il 



B' 1 Xi 

a x 2 



/ -1 2X' 



1 e 

4r? 2r r 



= 1 



(27) 



Fig. 5 and Fig. 6 are the graphs for R and M respectively for GCCG collapse. They almost resemble with 
Fig. 3 and Fig. 4 respectively. Only the 7? curve is almost a straight line whereas in case of MCG, it was a curved 
one. 



3.4 Effect of a combination of dark matter and dark energy(in the form of MCG) 
3.4.1 Case I : Q = i.e., No Interaction Between Dark Matter And Dark Energy: 

Energy density of matter as given by equation (|T4")) is pm — §r ■ Also the energy density of MCG is given in 
equation (|20l) . The expressions for relevant physical quantities are 



R(T) = -r s a 







I Co 


A 




a 3 



x 1 



Ci_ 

n x 2 



2X 3 - 



1 e 
472 + 27: 



and 



M(T) = -r|a 3 







\C 


A 


w 


a 3 



AA 



Ci_ 



2X 3 ' 



1 e 

4r 2 2r r 



The limiting values of the physical parameters are as follows, 

easel 

p M 00, Pmcg ->• 00, for 1 + A > 0; Pa/cg 



B 

TTa 



(28) 



(29) 



/or 1 + A < 
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Fig. 7 Fig. 8 

Fig 7 : The time derivative of the radius is plotted against time, r = 10, r c = 1000, Co = 
0.00001, d = 0.05, Xi = 40, X 2 = 0.04, X 3 = 0.4 is considered. 

Fig 8 : The mass of the collapsing cloud is plotted against time, r = 10, r c = 1000, Co = 0.00001, C\ = 
0.05, Xi = 40, X 2 = 0.04, X 3 = 0.4 is considered. 



Collapse with DE and DM in DGP(non-zero Q) Collapse with DE and DM in DGP(non-zero Q) 
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Fig 9 : The time derivative of the radius is plotted against time, r = 10, r c = 1000, n = 1, A = 
1, B = 1000, a = 0.5, zq = 1 is considered. 

Fig 10 : The mass of the collapsing cloud is plotted against time, r = 10, r c = 1000, n = 1, A = 
1, B = 1000, a = 0.5, zq = 1 is considered. 
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R 



—ear 
2r r ' 



for 1 + A > 0; M (T) 



i+a „-3A^,3 



C7 x l+a a 



-, /or 1 + A>0 



a ->• oo : pm -» 0, pa/cg — 



1 + A 



i? 



-v / 2ra-( 1+3A )C 1 2<1+Q) 



case2 

, /or 1 + A>0; pMCG^oo, for 1 + A < 



/or 1 + A<0; M{T) 



6 



/or 1 + A < 



3.4.2 Case II : Q 7^ i.e., Interaction Between Dark Matter And Dark Energy: 

In this case walking in the path of Cai and Wang [27j [28] we will assume the relation 



Pmcg 
Pm 



Ca 



3n 



(30) 



where C > and n are arbitrary constants. Now we solve the conservation equations §5§ and ^ to get the 
following expression for pr where pr = Pmcg + Pm, 



Pt +1 



(a + 1)5 (Ca 



[a (n - 1) - 1] fCa 3n + i\U a + 1 - na ')+ A ( a + 1 ") 
1 + a — not l + n + a + A + Act 1 + n + a — na Ca 3n 



Ca 3n (l + Ca 



3n \ 



n 

-(a+l) 



' 1 + Ca 3r 



(31) 



where 2-P1 is the hypergeometric function and zq is the integration constant. Using the relation px — Pmcg+Pm 
we get 

Ca 3n p T 



Pmcg 



1 + Ca 3n PM 1 + Ca 3n 

The expression for interaction in this case, is obtained by using equations HI [51 and 1121 and is given by, 



PT 



(32) 



Q = 



3pj 



'PT 1 e 
3 4r? 2r f 



1 + Ca 3 ' 

In this case the value for the gradient of scale factor is given by 



1 - l Ca?n (l + a 3n (C(l + A + 3n)) - BCT"^ 3 "" (l + 4 



3n\ Q + 1 -a-1 

" ) Pt 



(33) 



^3 

^"(l + Ca 3 ")- 4 



where A3 = 57 



3C 



(a + l)B 



a(n-l)-l 



(l + C*a 3 ")^( Q+1 -" Q ) 

C*a 3n 
1 + Ca 3n 

, and zi = zq 



2F1 



1 + a — na l + ?i + a + A + Aa 1 + n + a — na 





1 \ 


+ Z1 


+ 4r c 2 J 



e 
2^ 



(34) 



(a+l)-B 



The corresponding expressions for i? and 



mass M(T) is given as follows: 



1? 



-r E a 



-43 



a 3 »(l + Ca 3n ) A 



^^(a+l-na) fl 6(a+l-na) 

(l + Ca 3 ")* (Q+1 ~™ a) 



1 + a — na l+n + a+ A + Aa 1 + n + a — na 
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Ca 



3 n 



1 + Ca 3 ™ 





=** 1 \ 


+ 21 


+ 4r c 2 J 



e 

2rv 



(35) 



and 



M(T) = -a 3 r| 



2 -Pi 



a 3 ™(l + Ca 3 ") A 



^^(a+l-na) ^6(a+l-na] 

(l + Ca 3 «)*( Q+1 -" Q ) 



+ a — na 1 


+ n + 


a + 


A + Aa 1 + n -f 


- a — na 


n 




n 




n 


Ca 3 ™ 




i 


+— V+— 


2 


1 + Ca 3 ™_ 


+ 21 




4r;? J 2r c 





(36) 



From the above expressions the limiting behaviour of the physical parameters are obtained as follows: 



easel 



a —> : pm —> a 3 ™ , Pmcg a constant, pT — >• a Jn , R —> 



3n 



l_3n 

a 2 



M (T) — > a 3(1 - n) 



a — > oo 



Pm 



-3n(a+l)(A+l) 



Pmcg a 
M(T) 



case2 

3n(a+l)(A+l) 



— > a 



-3n(a+l)(A+l) 



,1-^(1+^) 



3[l-n(l+,4)] 



Fig. 7,8,9 and 10 are for the cases where DM and DE are present together. First two characterizes the no 
interaction scenario, whereas the last two deals with the interaction cases. 

When there is no interaction the magnitude of R decreases at first followed by a steep increment. This can 
be interpreted as the domination of DE over DM at the first phase and domination of DM over DE at the latter 
stage. In either case the mass decreases for no interaction case and lastly becomes asymptotic with time axis. 

For interacting DM and DE model, it is clear that DM dominates and the basic nature of the graphs are 
similar to DM collapse case. 

3.5 Effect of combination of dark matter and dark energy(in the form of GCCG) 
3.5.1 Case I : Q = i.e., No Interaction Between Dark Matter And Dark Energy: 



The energy density of DM is given by equation (jT4"]) and the energy density of GCCG is given by equation 
Using these two equations in the modified Friedmann equation for DGP brane, i.e., eauation ((T2|) . we obtain the 
expressions for the relevant parameters as given below. 

The time derivative of the geometrical radius of the collapsing star is given by, 



R(T) 



-dry, 



\ 




(37) 



The expression for mass of the collapsing cloud is given by, 
M(T) 



1 



a 3 r% 



\ 




(38) 



The limiting values of the physical parameters are given below, 

easel 
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Collapse with DE(GCCG) and DM in DGP(Q=0) 
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Fig 11 : The time derivative of the radius is plotted against time, r = 10, r c = 1000, n = 1, 
10, C = 0.00001, B' = 5, Xi = 0.00005, X 2 = 1.5, X 3 = 0.25 is considered. 
Fig 12 : The mass of the collapsing cloud is plotted against time, r = 10, r c = 1000, n = 1, 
10, C = 0.00001, B' = 5, X x = 0.00005, X 2 = 1.5, X 3 = 0.25 is considered. 
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Fig. 13 Fig. 14 

Fig 13 : The time derivative of the radius is plotted against time, r = 10, r c = 100000, n = 1, 
0.001, B' = 5, a = 0.5, C = 1000 is considered. 

Fig 14 : The mass of the collapsing cloud is plotted against time, r = 10, r c = 100000, n = 1, 
0.001, B' = 5, a = 0.5, C = 1000 is considered. 
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3.5.2 Case II : Q ^ i.e., Interaction Between Dark Matter And Dark Energy: 

Like the case of MCG here also we assume the relation 



Pgccg 
Pm 



= ka 3 



(39) 



where k > and n are arbitrary constants. Using this equation and equation [5U in the equations [5] and [Bl we 
calculate the expression for total energy density (pr) as given below, 
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(40) 



where C is the integration constant. The above relation is obtained by providing particular values to some of 
the variables involved in the equation, since the differential equation obtained was not solvable directly. Here 
the particular values used are lj = — h, k = 1, n = 1 and a = 1. Using the relation px = Pm + Pgccg, we get, 
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1 + ka 3n 

The expression for the time derivative of scale factor is given by, 
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The corresponding expressions for other physical parameters are, 
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(43) 



and 



M(T) = -a 3 r 3 
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(44) 



Using the conservation equations [5] and El and using the expressions for pm, Pgccg and Pt we get the 
following expression for the interaction, 




L-=- ( ka 3n (n + 1) + 1 - fc~ a f P f = 

l + ka 3n \ ' \l + ka 3n 
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ka in p T 



1 + fca 3r \ 

Fig. 11, 12, 13 and 14 are for GCCG case. These are quite similar with figs. 7, 8, 9 and 10 respectively. 




1 



(45) 



4 Gravitational Collapse in Loop Quantum Cosmology 

In recent years, loop Quantum Gravity (LQG) is an outstanding effort to describe the quantum effect of our 
universe |29l 130) . LQG is a theory trying to quantize the gravity with a non-perturbative and background 
independent method. The theory and principles of LQG when applied in the cosmological framework creates 
a new theoretical framework of Loop Quantum Cosmology (LQC) [3TJ [32j [33]. In this theory, classical space- 
time continuum is replaced by a discrete quantum geometry. The effect of LQG can be described by the 
modification of Friedmann equation to add a term quadratic in density. In LQC, the non-perturbative effects 

2 

lead to correction term — to the standard Friedmann equation. With the inclusion of this term, the universe 
pi 1 ' 

bounces quantum mechanically as the matter energy density reaches the level of p\ (order of Plank density). 

Recently the model of DE has been explored in the framework of LQC. The cosmological evolution in LQC 
has been widely studied for quintessence and phantom DE models 34 . The modified Friedmann equation for 
Loop Quantum Cosmology is given by ([351 IM1 157] ) 



Pt 



H = — \ 1 — — 

Pi 



Pt 



(46) 



Here p\ = ^/ii: 2 ^ 3 G 2 h is the critical loop quantum density and 7 is the dimensionless Barbero-Immirzi param- 
eter, pt = Pm + Pe represents the total cosmic energy density, which is a sum of energy density of DM (pm) 
and the energy density of DE (pe)- 

Like the previous section in this section also we will study the role of DM and DE separately during the 
collapse. Then we will study the effect of the combination of DE and DM in the collapsing process, both in the 
presence and in the absence of interaction. 
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Collapse with DM in LQC Collapse with DM in LQC 




Fig. 15 Fig. 16 



Fig 15 : The time derivative of the radius is plotted against time. i?o = 100, uq = 0.4, Nq = 0.1, T 
10 is considered. 

Fig 16 : The mass of the collapsing cloud is plotted against time. Rq = 100, no = 0.4, No = 0.1, T 
10 is considered. 



4.1 Collapse with Dark Matter 

Here pm ^ and pe — 0. From the conservation equations ([5]) and © we get the solution in (TT4")) . Putting the 
obtained value of Pm in equation (|46[) we get 



a 3 = [n Q (T - T)f + TV 2 



(47) 



where uq = y | Co and Nq = y J^J- and To is the integration constant. The time derivative of the geometrical 
radius of the collapsing object is determined as follows 



2 (a 3 - Na 2 ) 2 
R(T) = --R ± 

The total mass of the collapsing cloud as determined is given by 



(48) 



M(T) 



9 



-R, 



nod 3 



(49) 



where Rq = Hq. The limiting values of the physical parameters are given as T — >■ oo, a — >• oo, pa/ 



0, E -> 0, M (T) 



2J?.p 
9no ' 



Like the previous case, the time for formation of apparent horizon is given by 



the real root of the equation R 2 (Tah) = ^|a 2 (T4ff) = 1 as given by equation (fTTj). Thus the corresponding 
expression for LQC is, 



^Rlj [n (To - T)] 6 - [n 2 (T - T) 2 + TV 2 ] * 
From the above expression it can be seen that singularity will form at the instant T = Tq. 

4.2 Collapse with Dark Energy in the form of Modified Chaplygin gas 



(50) 



In this case, pu = 0. The expression of the density of MCG as given in equation (|20p is pmcg = j^x t a 3(i+A)(i+ a ) 
Using this relation in the modified Friedman equation for Loop Quantum Cosmology given in equation (|46[) , we 
get the expressions for the other relevant physical quantities as follows, 
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Fig 17 : The time derivative of the radius is plotted against time, r = 0.5, p\ = 1000, C\ = 0.5, X\ 
0.4, X 2 = 0.4, X 3 = 0.4 is considered. 

Fig 18 : The mass of the collapsing cloud is plotted against time, r = 0.5, p\ = 1000, C\ = 0.5, X\ 
0.4, X 2 = 0.4, X 3 = 0.4 is considered. 
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The total mass of the collapsing cloud as determined is given by 

x 3 



M(T) 



G 



X 4- Cl 

a X2 



Pi I a* 2 
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(51) 



(52) 



We know that for the cloud to undergo collapse, we should have R < 0, and hence from the above expression 
for i?,we have 



a > - 

for X2 > 0, i.e., 1 + A > 0. Again we should have, 

1 

a < 



(53) 



(54) 



for X2 < 0, i.e., 1 + A<0, where X2 = —p. The limiting values for the physical parameters are given as 
follows: 

easel 
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/3pi 
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Fig 19 : The time derivative of the radius is plotted against time, r = 1000, p\ = 1000, B' = 5, C 
100, X 1 = 0.00005, X 2 = 1.5, X 3 = 0.25 is considered. 

Fig 20 : The mass of the collapsing cloud is plotted against time, r = 1000, p\ = 1000, B' = 5, C 
100, X 1 = 0.00005, X 2 = 1.5, X 3 = 0.25 is considered. 
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The cloud will start to be trapped at the instant given by the real roots of the equation, 
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(55) 



4.3 Collapse with Dark Energy in the form of Generalised Cosmic Chaplygin Gas 

Here pm = 0. The expression for the density of DE is given by equation (|24[) . The expressions for the relevant 
physical quantities are, 



and 
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For the cloud to undergo collapse, we should have, 

. 1 



a > — 

Pi 



2X' 



(56) 



(57) 



(58) 
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Collapse with combination of DE and DM in LQC(Q=0) 



Collapse with combination of DE and DM in LQC(Q=0) 





Fig. 21 



Fig. 22 



Fig 21 : The time derivative of the radius is plotted against time, r = 10, p\ = 1000, B' = 5, C\ 
500, C = 0.00001, Xi = 40, X 2 = 0.04, X 3 = 0.4 is considered. 

Fig 22 : The mass of the collapsing cloud is plotted against time, r = 10, p\ = 1000, B' = 5, C\ 
500, C = 0.00001, X x = 40, X 2 = 0.04, X 3 = 0.4 is considered. 

The limiting values of the corresponding physical parameters are given as follows, 

easel 



a -> : Pgccg -> oo, for 1 + w > 0; Pgccg -> [C + 1]* 
R(T) 



for 1 + cj < 



r^a" 



nip: 



; for l+w>0 



M(T) 



fl (i+o0d+") , /or 1 + ui > 



case2 

a -> oo : pgccg oo, for l+w<0; pgccg -> [C + 1] 1 



i?(T) 



M(T) 



£>J~Pi 



(C + l)5tlT=j pi - (C" + l)i+° 



(C' + l)~ pi-(C' + l)~ 



/or 1 + cj > 
/or 1 + cj > 



for 1 + w > 



All the above limiting values are calculated considering that equation (|58| is satisfied. The time for the formation 
of the apparent horizon is given by the real root of the following equation, 
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(59) 



4.4 Effect of combination of Dark Energy(in the form of MCG) and Dark Matter 

4.4.1 Case I : Q = i.e., No Interaction Between Dark Matter And Dark Energy: 

The solution for pm and pmcg are given in equations (|14[) and (|20|) . The relevant physical quantities are 
obtained as 
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(60) 
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Collapse with DE and DM in LQC(non-zero Q) 
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Collapse with DE and DM in LQC(non-zero Q) 
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Fig 23 : The time derivative of the radius is plotted against time, r = 10, p\ = 10, B = 5, C 
5, , n = 1 is considered. 

Fig 24 : The mass of the collapsing cloud is plotted against time, r = 10, pi = 10, B = 5, C 
5, , n = 1 is considered. 
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In this case, the condition for the cloud to undergo collapse is given by, 



(61) 



(62) 



for X2 > 0, i.e., 1 + A > 0, and a — > 00. The above result is also valid for X2 < 0, i.e., 1 + A < 0, and 
a — > 0. The limiting values of the physical parameters are as follows, 
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4.4.2 Case II : Q ^ i.e., Interaction Between Dark Matter And Dark Energy: 

As in the case of DGP brane here also we will assume the relation pMGG — Ca 3n as given in equation ([50]) . The 
expressions for p^j, Pmcg an d Pt ar e given in equations (fH|) . (IM1) and (f5T|) respectively. The expression for 
interaction is given by 
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1 Pt Pt 
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1 - —) U 7^ ( 1 + « 3n (CQ. + A + 3n)) - BC-° a - 3m (l + ^a 3n ) Q+1 p^- 1 ) - 1 

/°i / Li + ^ a ^ ' 

(63) 
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In this case the value for the gradient of scale factor is given by 
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where pr is given by (|3"Tj) . The corresponding expressions for mass M(T) and R are given as follows: 
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The limiting values for the physical parameters for this case are given by 
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4.5 Effect of combination of Dark Energy(in the form of GCCG) and Dark Matter 
4.5.1 Case I : Q = i.e., Non-interaction Between Dark Matter And Dark Energy: 

The solution for pm and pgccg are given in equations (1141) and (|24p respectively. The relevant physical 
quantities are obtained as 
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(67) 



(68) 



The collapsing process will take place provided R{T) < 0, hence the corresponding condition is given by, 
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(69) 
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Fig. 25 Fig. 26 

Fig 25 : The time derivative of the radius is plotted against time, r = 100, p\ = 1000, B' = 
5, C = 100, , n = 1, Xi = 0.00005, X 2 = 1.5, X 3 = 0.25, C = 0.00001 is considered. 
Fig 26 : The mass of the collapsing cloud is plotted against time, r = 100, p\ = 1000, B' = 5, C = 
100, , n = 1, Xi = 0.00005, X 2 = 1.5, X 3 = 0.25, C = 0.00001 is considered. 



Collapse with DE(GCCG) and DM in LQC(non-zero Q) Collapse with DE(GCCG) and DM in LQC(non-zero Q) 




Fig. 27 Fig. 28 

Fig 27 : The time derivative of the radius is plotted against time, r = 10, p\ = 100000, B' = 5, C = 
0.001, n = 1 is considered. 

Fig 28 : The mass of the collapsing cloud is plotted against time, r = 10, p\ = 100000, B' = 5, C = 
0.001, n = 1 is considered. 
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The corresponding limiting values of different parameters are given below, 
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4.5.2 Case II : Q ^ i.e., Interaction Between Dark Matter And Dark Energy: 

As usual we consider PG p° G = ka 3n , where k > and 'n' is an arbitrary constant. Just like the previous cases 
the expressions for pu, Pgccg an d Pt are also taken into account from equations (TT^j) , (IM1) and P0|) 
respectively. The relevant parameters, i.e, d, i?(T) and M(T) are given by the equations (l64l) . (1651) and (|6G|) 
respectively using the values of pm, Pgccg and pt from the equations (ITU) , (jMf and respectively. 
Using the conservation equations (5) and (6), we calculate the interaction as given below, 
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(70) 



From figs. 15 to 28 we have demonstrated the R and mass curves for LQC. Here the R for DM collapse is 
increasing in magnitude with time. But the mass is decreasing. The fact that the velocity of inward collapse is 
increasing but mass is decreasing, is quite an awkward incident to explain. This is a very unique property for 
LQC. In figs. 17 and 18, we see the DE collapse in the universe filled up with MCG. The inward velocity with 
time whereas the mass of the collapsing centre decreases. The same result is followed for GCCG case in graphs 
figs. 19 and 20. Unlike DGP case in LQC DM and DE combined case is dominated by the speedy collapse and 
decreasing mass incident. 



5 Detailed Graphical Analysis 

In the above figures the time derivative of the geometrical radius(-R) of the collapsing object and the mass(M) 
of the same have been plotted against time, for all the above cases. From the plots it is evident that the time 
derivative of the radius inevitably takes negative values, which in turn shows that the radius decreases with 
time, thus resulting in gravitational collapse. We discuss the graphs for the two sections separately in detail. 

DGP BRANE MODEL 

From the first figure we see that there is a considerable decrease in the value of R, towards more negative 
values in the presence of dark matter(dust) only. This shows that matter field collapses due to their inward 
gravitational pull, and gradually increases in mass by accumulating more matter, which is evident from fig. 2. 
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From fig. 3, it is seen that in the presence of dark energy in the form of MCG, value of R increases with time 
and has a gradual tendency of stepping into the positive region. This is due to the fact that the large negative 
pressure of the DE does not facilitate collapse under normal conditions. Initially if the space-time is trapped, 
during evolution it gets un-trapped and eventually the cloud shows a tendency to expand due to the presence 
of DE in the form of MCG. Even if the cloud does not expand, due to its local condensation, the time for 
formation of the apparent horizon and consequently a BH will be much delayed in future. Hence the formation 
of BH is quite an uncertainty, if at all it is a possibility. Fig. 4 shows that due to DE accretion on the collapsing 
object, there is a gradual decrease in its mass, with time. With GCCG as the DE the story is almost the same, 
as is evident from figures 5 and 6. Fig. 7 shows that in the absence of any interaction, a combination of MCG 
and dust undergoes a sudden collapse, forming a singularity of the space-time. As expected there is a decrease 
in mass due to presence of DE in fig. 8. In the presence of interaction neither DM nor DE(MCG) can dominate 
over each other, and hence we observe a mean behaviour between the two in fig. 9. From fig. 10 it is evident that 
there is an increase in mass with time. The reason being that interaction between DE and DM alleviates the 
DE accretion phenomenon on the collapsing object, thus preventing the mass from decreasing with time. With 
GCCG as DE we get almost identical scenario as MCG, in the absence of interaction with DM, which is evident 
from figs. 11 and 12. In the presence of interaction, we get a slightly different situation in case of GCCG as far 
as mass is concerned. In fig. 14 we see that there is a gradual decrease of mass, thus showing the dominance 
of DE over DM, unlike the case of MCG. This is characteristic of GCCG. The above results are obtained by 
considering large values for r c , the cross over scale, since for small values of r c , there is no theoretical possibility 
for a collapse. 

LOOP QUANTUM COSMOLOGY MODEL 

From fig. 15 we see that there is at first a gradual decrease and then a steep decrease of mass due to collapse 
of DM in LQC model. This result is totally different from the corresponding result for DGP brane model. This 
is indeed quite surprising. The gravitational pull in case of LQC can be speculated to be really weak compared 
to other gravity theories. Fig. 16 is quite obvious. Collapse with MCG in LQC is exhibited in figs. 17 and 18. 
The results obtained tally with the results obtained for DGP brane model. The only difference being that in 
case of LQC the graph for R rises more steeply than in case of DGP brane. This again reflects on the weakness 
of gravity, speculated in the previous figs. Figs. 19 and 20 characterizes the collapse of GCCG in LQC. Here 
also the results are similar to that of DGP brane model, the only difference being the steepness of the slope 
in LQC just like the previous graphs. Figs. 21 and 22 shows the plots for combination of MCG and DM in the 
absence of interaction. Unlike DGP brane model in fig. 21 there is no sudden collapse. In fact the increase in 
value of R shows the dominating nature of DE over DM. The large negative pressure of MCG prevails over 
the weak gravity of DM, and hence there is a comparatively lesser tendency of collapse. Fig. 22 shows decrease 
of mass with time, due to MCG accretion on the collapsing cloud. In presence of interaction MCG and DM 
show much moderate behaviour as is evident from the figs. 23 and 24. Association with each other suppress 
their individual behaviour upto a great extent, thus bringing their mean properties to light. From figs. 25 and 
26 we get the characteristics of collapse with DM and GCCG as DE. Here also unlike DGP brane model there 
is no possibility of a sudden collapse, thus showing the dominance of DE over DM. Finally in figs. 27 and 28 we 
consider collapse with GCCG and DM in presence of interaction. Unlike the case of MCG, in this case there is 
a greater possibility of collapse as seen from the figs. 

Hence in case of LQC, the effect of MCG is felt upto a much greater extent than the effect of GCCG when 
combined with DM. This is characteristic of the two fluids and solely depend on their respective pressures. The 
above results are obtained on considering large values for p\ , the critical loop quantum density, since for small 
values of p\ there is no theoretical possibility of collapse. 

6 Discussions and Concluding Remarks 

Here we have studied the gravitational collapse of a spherically symmetric dust cloud of finite radius, filled 
with homogeneous and isotropic fluid. The study was carried out separately in two different types of modified 
gravity theories, namely DGP brane model and Loop quantum cosmology model in the background of a unified 
model of DE and DM. Two candidates for DE has been considered separately, namely modified Chaplygin gas 
and Generalized Cosmic Chaplygin gas. At first the effect of DM and DE(both MCG and GCCG) was studied 
separately and then their combined effect was studied both in the presence and absence of interaction. Values 
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for relevant parameters (time derivative of radius and mass of the collapsing cloud) were found out and their 
variations with time was plotted. A detailed graphical analysis was done to get a comparison between the 
results obtained in the two different models. 

It was seen that in case of DM, there is a high possibility for the matter cloud to undergo collapse, and 
the consequent formation of a BH was inevitable, for both DGP and LQC models. In the background of DE, 
i.e., MCG or GCCG, the cloud shows a reluctance to undergo collapse, due to the high negative pressure of 
the DE fluids. If at all there is any possibility of a collapse and subsequent formation of a BH, it will be 
much delayed in future. Hence the formation of BH is absolutely an uncertain phenomenon in this case. This 
phenomenon is much more pronounced in LQC model than in the DGP brane model, thus showing the relative 
weakness of gravity in the LQC model compared to the DGP brane model. This is a very interesting result 
indeed. In the unified model of DE and MCG in the absence of interaction, there is a sudden collapse showing 
the the dominance of matter over DE in the DGP model. Hence the fate of the collapse is a BH. But we have 
a completely different scenario for LQC model, where there is no such collapse showing the dominance of DE 
and the weakness of gravity. In the presence of interaction we see a moderation as far as the individual effects 
of DM and DE are concerned. The presence of either of them suppress the extremeness of their behaviour, as 
individuals. In all the above cases the presence of DE decreases the mass of the cloud due to DE accretion on 
it. In case of DM only there is an increase in mass due to accumulation of more matter. 
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